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Cu ^ The Hamilton-Jacobi problem is revisited bearing in mind the consequences arising from a 

possible bi-Hamiltonian structure. The problem is formulated on the tangent bundle for La- 
grangian systems in order to avoid the bias of the existence of a natural symplectic structure 
on the cotangent bundle. First it is developed for systems described by regular Lagrangians 
and then extended to systems described by singular Lagrangians with no secondary con- 
straints. We also consider the example of the free relativistic particle, the rigid body and 
the electron-monopole system. 
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1 Introduction 

Hamilton-Jacobi theory provides important physical examples of the deep connection between 
first-order partial differential equations and systems of first-order ordinary differential equations. 
In this respect, it is also a stepping stone to the Schrodinger wave equation in quantum mechan- 
ics, and takes us as close as possible, within classical theory, to the notions of wave function 
and state in quantum theory. As a matter of fact, we obtain Hamilton- J acobi- type equations 
whenever we consider a short-wave approximation for the solutions of wave- type equations, 
i.e., hyperbolic- type equations (this includes the classical limit of quantum mechanics in the 
Schrodinger picture by means of eikonal coordinates and the geometrical optics limit of wave 
optics) ^7]. Within the framework of wave mechanics, a complete solution of the Hamilton- 
Jacobi equation allows us to reconstruct an approximate solution of the Schrodinger equation 
by providing us with the phase of the wave function and the amplitude via the van Vleck deter- 
minant constructed out of the Hessian of the complete solution itself (see, for instance, ^Tj; P- 
172). 

As regards the Hamiltonian formulation of geometrical optics, one may recall that the origin 
of the whole method of canonical transformations in analytical mechanics can be traced back 
to the famous memoirs on Optics presented by Hamilton to the Royal Irish Academy. There, 
Hamilton showed that the propagation of wavefronts can be entirely characterized by the knowl- 
edge of a single function called the characteristic function. He also showed that the characteristic 
function obeys a first order partial differential equation, the so-called eikonal equation which is 
strictly related to the Hamilton-Jacobi equation. 

The name Hamilton-Jacobi is justified by the contribution given by Jacobi that the dynam- 
ical problem (the ordinary differential equation) is completely solved once a complete solution 
of the associated partial differential equation is known. 

Taking into account the fact that these equations were discovered almost two centuries 
ago, one may believe that everything must be known for them. As a matter of fact, we will 
argue and show in this paper, and in the forthcoming ones, that there are several aspects 
which have so far not been considered. Our own interest in reconsidering the Hamilton-Jacobi 
theory was generated by the existence of bi-Hamiltonian descriptions for completely integrable 
dynamical systems and the desire to unveil and understand the quantum counterpart of bi- 
Hamiltonian systems. In particular, due to the relevant role of the Hamilton-Jacobi theory in 
the Schrodinger picture, it seems appropriate to achieve a proper understanding of the Hamilton- 
Jacobi formulation for bi-Hamiltonian systems as a preliminary step toward the the possibility 
of a better understanding of the corresponding quantum situation. 

Vinogradov [H7j has exhibited a deep relation between the commutation relations of differ- 
ential operators acting on functions over the configuration space Q and the canonical Poisson 
brackets of their principal symbols on the cotangent bundle T*Q. This connection seems to rule 
out the possibility of considering the Hamilton-Jacobi version of the bi-Hamiltonian systems. To 
escape this apparent impossibility we find convenient to formulate the Hamilton-Jacobi theory 
on the tangent bundle TQ with the help of a regular Lagrangian function and the associated 
Lagrangian two-form. Thus we remove the bias of a natural symplectic structure on our carrier 
space, unlike in the case of the cotangent bundle. Working with Lagrangians on TQ we have 
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the possibility of dealing more directly with relativistic aspects and with dynamical systems 
described by degenerate Lagrangians (gauge theories), and therefore with the classical limit of 
their corresponding quantum systems. 

In this paper we will not address the problem of constrained Lagrangian dynamics in full 
generality; this will be done in a forthcoming paper. Subsequently we also shall consider the 
classical limit of quantum bi-Hamiltonian systems and extend our Hamilton-Jacobi picture to 
classical field theories which allow for bi-Hamiltonian descriptions. 

As a spin-off from our tangent bundle formulation of the Hamilton-Jacobi problem, we will 
identify two main geometric aspects of the classical formulation: the first one consists of finding 
a foliation transverse to the fibers (of TQ or T*Q) and invariant under the dynamical evolution, 
while the second one requires that the foliation be Lagrangian with respect to a dynamically 
preserved symplectic structure. In this approach the dynamics (ordinary differential equation) 
plays a prominent role because we consider alternative Lagrangian or Hamiltonian descriptions. 
Therefore, our generalization is to search for invariant foliations of the carrier space with leaves 
having the same dimension as the configuration space Q, since we drop the requirement of 
"Lagrangianity" . Thus the partial differential equation associated with our problem (equation 
© in Section 121) will be a partial differential equation for a vector valued function rather than for 
a scalar valued function as in the standard formulation. The transition from the vector valued 
function to the scalar valued one takes place with the help of the symplectic structure which 
allows us to associate a closed 1-form (and therefore locally a function) with our vector field by 
requiring the foliation to be Lagrangian. 

When considering geodetical motions on Lie groups, an interesting situation arises in which 
the first step is accomplished but the second one is problematic, as we will see in Section 7. 
Similar aspects also emerge when dealing with bi-Hamitonian systems or systems described by 
equivalent Lagrangians: Here we find invariant foliations; they may be Lagrangian with respect 
to one symplectic structure but not Lagrangian with respect to some other invariant symplectic 
structure. We shall discuss a few very simple examples to illustrate what is taking place. 

The paper is organized as follows. In Sections |2l and 01 we state the Lagrangian and Hamil- 
tonian geometrical formulations of the Hamilton-Jacobi problem, respectively, showing how the 
standard classical problem is a particular case of the extended one, and clarifying the geometri- 
cal meaning of particular and complete solutions |27| l,S5j . The relation between both formalisms 
is also discussed. Section |1] is devoted to extending the theory to the particular case of singular 
dynamical systems: those where there are no Lagrangian constraints or, what is equivalent, 
when secondary Hamiltonian constraints do not appear. As an application of the above case, 
the Hamilton-Jacobi problem for non-autonomous Lagrangian and Hamiltonian systems is dis- 
cussed in Section |S1 Finally, as examples, we apply our theory to the free relativistic particle in 
Sectional to the free motion on a Lie group, to the rigid body, and to the electron-monopole 
system in Section |3 

Notation: Throughout this paper Q is a n-dimensional differentiable manifold representing 
the configuration space of a dynamical system, and tq : TQ —> Q and ttq : T*Q -^ Q are 
its tangent and cotangent bundles, representing the phase spaces of velocities and momenta, 
respectively. 
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On the cotangent bundle there is a canonical symplectic form oj = —d9, where 6 is the 1- 
forni 6 with coordinate expression 9 = pidq^. Here and in the rest of this paper, sum over paired 
covariant and contravariant indices is understood. This symplectic form associates a vector field 
Zh to every function H £ C°°{T*Q), as the solution of the equation i{ZH)io = dH (see e.g. P 
for the details). 

On the tangent bundle, the canonical object is the vertical endomorphism S, with coordinate 
expression S = (d/dv^) (8> dq^. Given a Lagrangian function L £ C°°{TQ), we define the Cartan 
l-form 6l = dgL = {dL/dv^)dq^ and the Cartan 2-forni ul = —dOi- The dynamical vector 
fields associated to the Lagrangian are the solutions of the dynamical equation i{r)coL = dE^, 
where El = A(L) — L £ C°°{TQ) is the Lagrangian energy function and A G X{TQ) is the 
Liouville vector field (see e.g. j23l I12j and references therein). 

We also remark that T{TQ) has two different vector bundle structures over TQ, given 
respectively by ttq: T{TQ) -^ TQ, i.e. considering TQ as new configuration space, and 
Ttq : T[TQ) -^ TQ. Maps X : Q ^ TQ that are sections for tq are the vector fields in Q, and 
the set of such vector fields will be denoted by 3i{Q). Correspondingly, maps X: TQ — * T{TQ) 
that are sections for ttq are the vector fields in TQ, and those which are also sections for Ttq are 
said to be second order differential equation fields (hereafter referred to as SODE vector fields). 
This means that their integral curves, which are the trajectories of the system, are holonomic. 
A vector field X £ 3i{Q) can be lifted to TQ producing the so called complete or tangent lift of 
X and denoted by X^ £ X{TQ). More details can be found in |131l3Uj . 



2 Lagrangian formulation of the Hamilton— Jacobi problem 

Li this section we formulate the Hamilton-Jacobi problem on the tangent bundle. In this setting 
we are able to handle dynamical systems which admit alternative Lagrangian descriptions, and 
we clearly show how the search for solutions of the Hamilton-Jacobi problem splits in two steps. 

We recall that, in the standard formulation, the Hamilton-Jacobi problem consists in finding 
a function S{t,q), known as the principal function, such that the partial differential equation 
(pde) 

ds ^( ds\ ^ 

is satisfied. If we put S{t, q) = W{q) — tE, where £■ is a constant, then the function W, known 
as the characteristic function, has to satisfy 

.(,.f)^.. 

Both of the above pde are known as the Hamilton-Jacobi equation. However we will always 
refer to the second one. 

In geometric terms, equation (^ can be written as {dW)*H = E, where we understand dW 
as a section of the cotangent bundle. In other words, we look for a section a of T*Q such that 
a*H = E and a is a closed l-form, da = 0, and hence locally exact, a = dW. The second 
condition, da = 0, can alternatively be expressed in terms of the canonical symplectic form on 
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T*Q in the form a*ij = 0, so that one can reformulate the Hamilton-Jacobi equation in the 
form pQ 

a*H = E, a*uj = 0. (2) 

Consider now the Lagrangian formahsm. Let L G C^{TQ) be the Lagrangian function 
and 9l, ujl be the associated Cartan forms. A Uteral translation of the above coordinate- free 
formulation of the Hamilton-Jacobi equations from the cotangent bundle to the tangent bundle 
would be ini 

X*{El)=E, X*{ujl)=0, (3) 

where X : Q ^ TQ is the unknown "vector valued" function, and the second equation states 
that the vector field X is associated (at least locally) with a function W by means of the relation 
^*(^l) = dW, which is a stronger version of X*{ujl) = 0. 

Among the many important consequences that may be deduced from the existence of a 
solution of the Hamilton-Jacobi equation, let us recall the following. Let P{q) = ^^{q), and 
consider the vector field X = ^{q,P{q)). If q = j{t) is a solution of the differential equation 
q = X{q), then X{t) = {'y{t),P{'j{t))) is a solution of the Hamilton equations. The Lagrangian 
counterpart of this property reads as follows. If X is a solution of (j31) and q = 7(t) is a 
solution of the differential equation q = X{q), then ^{t) = ('j{t),X('j{t))) is a solution of the 
Euler-Lagrange equations. 

This fact will be our starting point in the study of the Hamilton-Jacobi equation and its 
generalization. We will look for the implications of this property and its relation with equa- 
tions (121). 

2.1 Statement of the problem and solutions 

We will assume first that the Lagrangian L is regular, and we will leave for Section [l] the 
analysis of the unconstrained singular case. The regularity of the Lagrangian is equivalent to 
the regularity of the Cartan 2-form, so that lol is symplectic. It follows that there exists a 
unique solution T^ G X{TQ) of the Lagrangian dynamical equation 

i{rL)uJL = dEL . (4) 



Tl is called the Lagrangian vector field of the Lagrangian system. It is well known |12j that F^ 
is a second order differential equation. 

Generalized Lagrangian Hamilton Jacobi problem. Let L G C°°{TQ) he a Lagrangian 
function. The generalized Lagrangian Hamilton-Jacobi problem consists in finding a vector field 
X : Q ^> TQ such that, if ^■. M. ^ Q is an integral curve of X , i.e. ^ = X o^, then 7: M ^ TQ 
is an integral curve ofT^; that is, 



X o 7 = 7 =^ Fi o 7 = X o 7 . 
X is said to be a solution of the generalized Lagrangian Hamilton-Jacobi problem. 

As we will see in a moment, in geometrical terms, this requirement means that the image of 
X, as a map from Q to TQ, is a F^-invariant submanifold of TQ. Let us show first an example. 
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Example 1 The dynamics of the free particle in M? is given by the regular Lagrangian function 



L{q\q^v\v') = l[{v'f + {v'f] , 



with associated geometrical objects 



9l = TT-r <iQ^ + TT^ dq^ = v^ dq^ + v^ dq^ 



2\2 
I 1 ■(/ I -1- I ■(/ 



The vector field 



El = ^[iv'f + (v , ^ 

UL = dq^ A dv^ + dq^ A dv^ 
id 2 ^ 
Qqi Qq2 



9 kq^ -I d 

X — K 7-— j- H ^ 7—^ , K,L €z 

oq^ q^ oq'^ 



defines a two-parameter family of vector fields on Q = 'M? which are generalized solutions. We 
also find that 

X*{u:l) = -'^^dq'Adq\ 



,2 fkq^ -I 
k^ + ' 



X*{El) = \ 



Thus, the simple translation of the geometrical relations from T*Q to TQ would be violated. 
Now we can formulate on TQ a PDE which replaces the pde for the characteristic function W . 
We find that it must be stated in terms of a vector valued function. 

Proposition \ X is a solution of the generalized Lagrangian Hamilton-Jacobi problem if, and 
only if, X and Tl are X-related; that is, 



TloX = TX oX . (5) 



Proof X is a solution of the generalized Lagrangian Hamilton-Jacobi problem if, for every 
7 : M ^ Q such that X o 7 = 7, then 



r^ 07 = ^07 = TX 07 = TX 0X07. 

But r^ o 7 = r^ o X o 7, and as X has integral curves through every point q & Q, this is 
equivalent to TX o X = Tl o X. 

The proof of the converse is straightforward. ■ 

This equation for a given SODE T^ defines a pde for X and replaces the pde for W in the 
standard formulation of the Hamilton-Jacobi problem. 

In addition we have: 
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Proposition 2 X is a solution of the generalized Lagrangian Hamilton- J acohi problem if, and 
only if, the submanifold Ini X C TQ is invariant by the Lagrangian vector field Tl (that is, Tl 
is tangent to the subm^anifold X{Q)). 

Proof For the direct implication, it suffices to show that, for every q £ Q, TL{Xg) is tangent 
to ImX, and it holds because, by proposition^ TL{Xg) = TqX{Xq). 

Conversely, if F^ leaves IraX invariant, then Ti^Xg) G Tx^^mX. Therefore, there exists 
u £ TqQ such that Fi(Xq) = TgX{u); hence 

Xq = {Ttq O TL){Xq) = {TqTQ O TqX){u) = Tq{TQ O X){u) = U , 

because tqoX = Hq, and F/, being a SODE, it is a section of the projection Ttq^ so Ttq oTl = 
Htq. Thus TL{Xq) = TqX{Xq) for every q G Q; that is, T^o X = TX o X, and X is a solution 
of the generalized Lagrangian Hamilton-J acobi problem by proposition ^ ■ 

If X is a solution of the generalized Lagrangian Hamilton-J acobi problem, then the integral 
curves of X are the TQ-projection of integral curves of F^ contained in ImX. 

Observe that we have not used that F^ is the Lagrangian vector field, so these results 
actually hold for every SODE F E X{TQ). Using the fact that F^ is the Lagrangian vector field 
of a Lagrangian system, the above results can be related with the energy Lagrangian function El 
in the following way, which avoids the explicit calculation of the dynamical Lagrangian vector 
field. 

Theorem 1 X is a solution of the generalized Lagrangian Hamilton-J acobi problem if, and only 

tf, 

i{X){X*uJL) = d{X*EL) (6) 



Proof From the Lagrangian dynamical equation Q we obtain 

X*i(TL)iOL = X*dEL = d{X*EL) , 
but, as X and Tl are X-related (proposition P) , we have that 

X*i(FiK = i(X)(X*cuL) , 

which yields ©• 

Conversely, suppose that X satisfies ©• The deviation Dl from the relatedness 

Dl = TloX-TXoX:Q^ ttq , 

is a vector field along X. We have to prove that Dl = 0. First we have that Di is rg-vertical. In 
fact, tqoX = Idg, and Tl being a SODE, it is a section of the projection Ttq, so TtqoTl = Utq, 
hence 

TtqoDl = TtqoTloX-TtqoTXoX = X-X = . 
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Furthermore, from the Lagrangian dynamical equation Q) we have X*i{TL)ujL = X*dEL = 
d{X*EL), which combined with the hypothesis, i(X)(X*a;L) = d{X*EL), leads to X*i{rL)ujL — 
i{X){X*ujL) = 0. Therefore, for every q £ Q and Yg G TgQ, we have 

= {X*i{TL)ojL-i{X){X*ojL))g{Y,) = iu;L)xi,){rLiq),TgXiYg)) - iX*ojL),{Xg,Y,) 
= {iOL)xi,){rL{q),T,X{Y,)) - {iOL)xi,)mX{Xg),TgX{Y,)) 
= iu;L)xi,){DLiq),T,X{Yg)) . 

Moreover, for every rg-vertical vector field V G X{TQ) we have that 

{u;LU,)iDLiq),ViX{q))) = -d(0i),(,)(Z)i(g), F(X(g))) , 

which vanishes for every q € Q. We recall that if a is a semibasic form, then da(Vi,V2) = 
for every pair of vertical fields Vi and V2. But 9l is a TQ-semibasic form, and Tx(q)TQ = 
Tx(q)iJrn.X) © Vx(g)(TQ), thus we have proved that 

{oJL)a{q){DL{q), Z{X{q))) = , for every g S Q, Z G X{TQ) , 

and hence DL{q) = 0, for every q £ Q, since ujl is nondegenerate. ■ 

To solve the generalized Lagrangian Hamilton-Jacobi problem is, in general, a hard task; it 
amounts to finding F/^-invariant submanifolds of TQ which are transverse to the fibers. Thus, 
it is convenient to consider a less general problem, which constitutes the standard version of the 
Lagrangian Hamilton-Jacobi problem: 

Lagrangian Hamilton— J acobi problem Given a Lagrangian function L G C°°{TQ), the 
Lagrangian Hamilton-Jacobi problem consists in finding solutions X of the generalized La- 
grangian Hamilton-Jacobi problem satisfying that X*ujl = 0. 

As = X*ujl = —X*d9L = —d{X*9L), we have that every point has an open neighborhood 
U C Q where there is a function W £ C°^(U) such that X*6l = dW (in U). 

Remark: In the example of the free particle in M'^ given above, the pull-back of the sym- 
plectic 2-form ujl by the vector field X is different from zero. Hence, X does not provide a 
solution of the Lagrangian Hamilton-Jacobi problem because it is not associated with a closed 
1-form on the configuration space. 

A straightforward consequence of the last theorem is: 
Corollary 1 If X is a solution of the Lagrangian Hamilton-Jacobi problem, then d{X*EL) = 0. 



Observe that if X is a solution of the Lagrangian Hamilton-Jacobi problem, then Im X is a 
Lagrangian submanifold of (TQ,ujl) contained in a level set of El- In fact, dim ImX = n and, 
if jx : ImX "^^ TQ denotes the natural embedding, we have that Jx^L = 0, due to X*{ujl) = 0. 

We can summarize the above results in the following: 
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Proposition 3 Let X £ 3l{Q) satisfy X*uji = 0. Then, the following assertions are equivalent: 

1. X is a solution of the Lagrangian Hamilton- J acohi problem. 

2. d{X*EL) = 

3. IraX is a Lagrangian suhmanifold ofTQ invariant by Tl. 

4- The integral curves ofT^ with initial conditions in IniX project onto the integral curves 
ofX. 



Coordinate expressions Let us show the local expressions of the objects so far presented. 
Consider coordinates (g*) on Q, and the corresponding natural coordinates (g*, w*) on its tangent 
bundle. 

Consider an arbitrary vector field T G X{TQ) satisfying the second-order condition, r(g, v) = 
{q,v;v,a{q,v)), and a vector field X £ X{Q): X{q) = {q,w{q)). Then we have 

iTXoX-To X){q) = (q, w{q); 0, ^w - a{q, w{q))\ , (7) 

which is a vertical vector field along X. Its vanishing is the necessary and sufficient condition 
for X and T to be X-related: 

-^w^{q)-a\q,w{q)) =0. 

This equation is the PDE for the vector valued function w^{q) which replaces the standard PDE 
for the scalar function W . 

dL d^L 
When r is the Lagrangian vector field F^, its components satisfy WijO^ = -—^ — . . v\ 

oq^ ov^ oqJ 

where Wij = . — r is the Hessian matrix of L. 
Then we can compute the 1-form 

dq\ (8) 



. ^, . Y*(.]p \ - I ^^^ ^^^ k 9"^^ j 9L 

\dv^dvi dq^ dv'^dqi dq^ 



v=w{q) 



whose vanishing also expresses that X is a solution of the generalized Lagrangian Hamilton- 
Jacobi problem. 

Looking carefully at the local expressions one can find a relation between —ixX*{uJL) + 
X*{dEL) and TX o X — T o X^ which is given by the Hessian, as we are going to show. 

To this end, let us first recall that, for any vector bundle E ^ Q, we have the vertical lift 
map vlg: E Xq E -^ YE C TE, an isomorphism which in fiber coordinates reads vl{q,u,v) = 
{q, u; 0, v). With E = TQ, this gives an isomorphism yI:TQxqTQ^ V{TQ) C T{TQ). 
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Associated with the Lagrangian L, we have the Legendre transformation 3'L: TQ -^ T*Q, 

dL 
which in coordinates reads 3'L{q,v) = (q^, ^-^). In other words, HL is the fibre derivative of L. 

Moreover, we can define the fiber Hessian S^L: TQ -^ T*Q (8) T*Q which defines, if the 
Lagrangian is regular, another isomorphism 3^L: TQ xq TQ -^ TQ xqT*Q. In coordinates, 
3^L{q,v) = {q\Wij) and ^L{q\v\u^) = {q\v\Wijui). 

With these ingredients, we achieve an alternative understanding of Theorem ^ 

Proposition 4 Let vl he the vertical lift map of the tangent bundle TQ. For any vector field X 
on Q, we have 

{X, -ixX*{ujL) + X*{dEL)) = ^o vpi o {TX oX -TloX). (9) 

Therefore, —ixX*{ujL) + X*{dEL) vanishes if and only ifTX o X — Tlo X vanishes too. ■ 

2.2 Complete solutions 

The most useful and essential idea in the standard Hamilton-Jacobi theory consists in finding, 
not only one particular solution as we have used in the previous subsection, but rather a complete 
solution of the problem. This may be defined as follows. 

Definition 1 Consider a solution Xx depending on n additional parameters A S A, where 
A C M" is an open set, and suppose that the map $: Q x A ^ TQ given by ^{q,X) = Xx{q) 
is a local diffeomorphism. In this case {Xx;\ G A} is said to be a complete solution of the 
generalized Lagrangian Hamilton-Jacobi problem. 

From the definition, it follows that a complete solution provides TQ with a foliation trans- 
verse to the fibers, and that the Lagrangian vector field T^ is tangent to the leaves. 

If {Xx]\ € A} is a complete solution, the integral curves of Xx^ for different A € A, will 
provide all the integral curves of the Lagrangian vector field Tl. This means that, if {qo,vo) G 
ImX, then there is Aq G A such that XxQ^qo) = vq, and the integral curve of Xx^ through qq, 
lifted by Xx^ to TQ, gives the integral curve of Tl through {qQ,VQ). This justifies the name of 
"complete solution" . 

Remark: We may use instead a fiber bundle P over A, such that F^ projects onto the null 
vector field; i.e. A is a space of constants of the motion and fibers have the same dimension as 
the configuration space Q. Thus we may take into account the nontriviallity of P as a bundle. 
On the other hand, if A were contractible the bundle would be trivial, and we would revert to 
the previous situation. 

Furthermore, different transversal foliations of TQ, with F^, tangent to the leaves, are dif- 
ferent ways to collect solutions of Ti smoothly and such that they project onto Q in a coherent 
way: integral curves of F^ in lva.Xx project onto integral curves of the associated vector field 
Xx. 

The relation between F^ and complete solutions is the following: 
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• If we have a family of n first integrals /i, ...,/„ of F^ such that dsfi A ... A dsfn 7^ 0, then 
ft = Ci, Ci £ R, for i = 1, . . . , n, define a transversal foliation. Thus we can locally isolate 
the velocities as functions of the coordinates g* and the constants q. Now, replacing in the 
expression of Tl these velocities and projecting to the basis, we obtain a local complete 
solution X(^, ,. \. 

• Conversely, if <I> : Q x A ^ TQ is a complete solution, then the functions defining locally 
the foliation give us the above family of integrals of motion of F^. More explicitly, the 
components of the map F: TQ -^ A given hy F = pr2 o <I>~^, are constants of the motion. 

Moreover, if the foliation is Lagrangian in (TQ,ujl), then we have a complete solution of 
the Lagrangian Hamilton-Jacobi problem. In this case the above family of first integrals are in 
involution. 

In our previous example of the free particle, varying the parameters {k, I) £ M? we obtain a 
complete solution. 

All these considerations are shown in the following example. 

Example 2 Let us consider the example of the two-dimensional standard harmonic oscillator 
described by 

L = liiv^ + {v^ - {qr - {e?) ■ 
The dynamical vector field is 

1^ 2 ^ 1 ^ 2 ^ 

Qql Qq2 Qyl Qy2 

and the standard Lagrangian symplectic 2- form is lol = dq^ A dv^ + dq'^ A dv'^. 
We know that the functions 

h = v\' + q\\ f2 = {v^f + (g^)^ h = (vY + iir, h = q'v' - q\' 

are constants of the motion. Of course, not all of them are functionally independent. Suppose 
their values are /i = C, /2 = 2Ei, f^ = 2E2, fi = L We can use, for instance, /2 and /a 
to express v^ and v'^ as functions of the base coordinates and the two parameters Ei and E2, 
and using these expressions in the dynamical vector field we find a vector field on the base Q 
depending on the two energies: 

Note that the two functions we have used are in involution, {/2, /s} = 0, and that 

{XE„E,yiOL = dq^ A d(V2^i-(gi)2) + dq^ A di^/2 E2 - {q^^) = . 

This is a 2-parameter family of vector fields, for which the images are Lagrangian submanifolds 
with respect to ujl- 
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We can also choose the functions /i and /4 for obtaining expressions of the velocities in 
terms of positions, when q^v"^ + q'^v^ 7^ 0, because 

dsfl A dsh = {v^ dq^ + v^ dq^) A {-q^ dq^ + q^ dq^) = {q^v^ + q'^v^) dq^ A dq"^ , 

and in this case, 



1 _ -l±y/p+Aq^q^{C-q^q'^) ^ _ I ± y^l"^ + 4q^ q^{C - q^ q^) _ I + q^v^ 



V = :: , V 



2g2 ' 2gi gi 

and we have the vector field in Q 



XcAq\q 



~ \ 2g2 1 9gl 



n±^P^4gig2((g_gig2y \ g 

However, notice that because of {/i,/4} = /2 — /s, we find (Xc^i)*a;i / 0. Therefore, Xc; 
is a complete solution for the generalized problem, but not for the standard Hamilton-Jacobi 
problem. 



3 Formulation of the Hamilton— Jacobi problem on T*Q 

3.1 Statement of the problem and solutions 

We now consider the Hamiltonian formalism in the cotangent bundle. Let H £ C'^{T*Q) be 
a Hamiltonian function, and denote hy to = —dO G ^'^{T*Q) the canonical symplectic form. 
There exists a unique vector field Zfj G X{T*Q) whose integral curves are the trajectories of 
the system; that is, the solutions of the Hamilton equation. Geometrically this means that Zh 
is the solution of the Hamiltonian dynamical equation 

i{ZH)i^ = dH . (10) 

Zf{ is called the Hamiltonian vector field of the system. 

As in the Lagrangian formalism, let us start with the generalized version of the Hamilton- 
Jacobi problem, which can be stated as follows: 

Generalized Hamiltonian Hamilton Jacobi problem Given a Hamiltonian vector field 
Zh G X(T*Q), the generalized Hamiltonian Hamilton-Jacobi problem consists in finding a 
vector field X : Q ^ TQ and a 1-form a: Q ^ T*Q such that, if ^■. W ^ Q is an integral curve 
of X, then a o j: M ^ T*Q is an integral curve of Zjj- That is, 

X o ^ = j =^ 007 = Zh o (a o 7) . (11) 



The first result is: 
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Proposition 5 Given a vector field Zh G X(T*Q), {X,a) satisfies the condition ill]) if, and 
only if, the vector fields X and Zu are a-related; that is, 

ZHoa = TaoX . (12) 



Proof If {X, a) satisfies the condition (|11() then, for every 7: M ^ Q such that X o 7 = 7, we 
have 

Zh 0007 = 007 = Ta 07 = Ta 0X07; 

but, as X has integral curves through every point q £ Q, this condition is equivalent to Zh^ol = 
TaoX. 

The proof of the converse is straightforward. ■ 

In fact, both elements {X, a) satisfying the condition (|TT|) are related, since, by composing 
both sides of the above equation (^21) with Tttq, and taking into account that ttq o a = Idg, we 
have the following immediate consequence: 

Corollary 2 // (X, a) satisfies the condition ill]) then 

X = Tttq o Zh o o- 



It is interesting to remark that we also have the following relation between X and a: 

X = JHoa , 
where S'H : T*Q — > TQ is the fiber derivative of the Hamiltonian function. 

In terms of our previous geometrical formulation, this amounts to X*{9l) = a when L is the 
Lagrangian function associated with H. 

As X is determined by a, we introduce the following: 

Definition 2 A solution of the generalized Hamiltonian Hamilton- J acohi problem, for Zh is a 
1-form a G il^(Q) such that, if ^i M. ^ Q is an integral curve of X = Tttq o Zh ° a, then 
Q o 7: M ^ T*Q is an integral curve of Zh; that is, 

Tttq o Zh 0007 = 7 =^ a o J = Zh o (a o 7) . 

Then X = Tttq o Zh o a is said to be the vector field associated with a. 

Example 1 (continued) 

Consider the Hamiltonian function for a free particle in M^ 

H{q^,q'^,Pi,P2) = - {pi^ +P2^) ■ 
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1 Id 

The 1-forni a = —r dq and its associated vector field X = — tA^itt-^t provide a solution for the 
q^ q^ oq'^ 

generalized Hamiltonian Hamilton-Jacobi problem, but we will see later that they do not give 
rise to any solution of the standard Hamilton-Jacobi problem. 

Proposition 6 Given a vector field Zh E X(T*Q), a 1-form a G Q,^{Q) is a solution of the 
generalized Hamiltonian Hamilton-Jacobi problem if, and only if, the submanifold Ima C T*Q 
is invariant under the flow of the vector field Zn (that is, Zn is tangent to the submanifold 
Ima). 

Proof If a G 0,^{Q) is a solution of the Hamiltonian Hamilton-Jacobi problem, and X = 
Tttq o Zh o a, then Zh o a = Ta o X, and thus Z//(a(g)) = Ta{X{q)), for every q £ Q. Hence, 
Zh is tangent to Ima. 

Conversely, if Ima is invariant by Zh then ZH{oL{q)) G T^iq^lma, which implies that there 
exists u G TqQ such that Zniaiq)) = Tga{u). Defining X by Tqa{Xq) = ZH{a{q)), then 
X is differentiable since X = Tttq o Zh o a. Hence X is a vector field in Q which satisfies 
Zh o a = Ta o X, and then a is a solution of the Hamiltonian Hamilton-Jacobi problem. ■ 

If a G il^{Q) is a solution of the generalized Hamiltonian Hamilton-Jacobi problem, taking 
into account Corollary [2 we can conclude that the vrg-projection of the integral curves of Zh 
contained in Ima are the integral curves of X. 

Observe also that until now we have not used that Zh is a Hamiltonian vector field, so these 
results actually hold for every vector field Z G 3i{T*Q). When Zh is the Hamiltonian vector 
field of a Hamiltonian system, the above results can be expressed in terms of the Hamiltonian 
function. 

As in the Lagrangian case, we can obtain an equation not involving directly the dynamical 
vector field: 

Theorem 2 Given the Hamiltonian vector field Zh G X{T*Q), a 1-form a G 0,^{Q) is a solution 
of the generalized Hamiltonian Hamilton-Jacobi problem if, and only if, 

i{X)da = -d{a*H) , (13) 

where X = Tttq o Zh oa is the vector field associated with a by means of the Hamiltonian vector 
field Zh- 

Proof From the Hamiltonian dynamical equation H1U() for Zh we obtain 

a*i{ZH)uJ = a*dH = d{a* H) . 
Furthermore, 6 is the canonical form of T*Q, so a*9 = a, and then 

a*uj = -a*de = -d{a*9) = -da , (14) 

therefore, as X and Zh are a-related, we have 

a*i{ZH)uJ = i{X)a*LL! = —i{X)da , 
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which yields (fT?^ . 

To prove the converse, first let us define 

DH = ZHoa-TaoX:Q^ TT*Q , 

which is a vector field along a. We have to prove that Dh = 0. First we have that Dh is 
TTQ-vertical; in fact, as ttq o a = Idg, 

Tttq o Dh = Tttq o (Z/f o a — Ta o X) = Tttq o (Z// o a — Ta o Tttq o Z// o a) 
= Tttq o Z// o a — Tttq o Z/f o a = . 



Furthermore, from the Hamiltonian dynamical equation (|lflj) and the hypothesis, as a*u> = —da, 
we have the following relations: 

a*i(Z//)a; = a*dH = d{a* H) , 
i{X)a*uj = -i{X)da = d{a*H) , 

and hence a*i{ZH)i^ — i{X)a*uj = 0. Therefore, for every q £ Q and Yq £ TqQ, we have 

= {a*i{ZH)uJ - i{X)a*Uj)q{Yq) = U)a{q){ZH{a{q)),Tqa{Yq)) - UJa(q){Tqa{Xq),Tqa{Yq)) 
= ^a{g){DH{q),Tga{Yq)) . 

Moreover, as V(7rQ) (the vrg-vertical subbundle oiTT*Q) is a Lagrangian distribution in {T*Q, cu), 
for every vrg-vertical vector field V £ X(T*Q) we have that 

'^a{q){DH{q),V{a{q))) =0 ; , 

for every q £ Q. But Ta(q)T*Q = r„(^)(Ima) © ^a(q){'^Q)i hence we have proved that 

^aiq){DH{q),Z{a{q))) =0 

for every q £ Q and Z £ X{T*Q); since a; is non-degenerate, we conclude that Dh = 0, or what 
is equivalent, X and Zh are a-related, and thus a is a solution of the generalized Hamiltonian 
Hamilton-Jacobi ■ 

As in the Lagrangian case, in general, to solve the generalized Hamiltonian Hamilton-Jacobi 
problem is a difficult task. So it is convenient to consider the following less general problem, 
which constitutes the standard version of the Hamiltonian Hamilton-Jacobi problem: 

Hamiltonian Hamilton Jacobi problem Given a vector field Zh £ X{T*Q), the Hamil- 
tonian Hamilton-Jacobi problem consists in finding a solution a £ Q,^{Q) of the generalized 
Hamiltonian Hamilton-Jacobi problem which is moreover closed, da = 0. 

As a consequence, every point has an open neighbourhood U C Q, where there is a function 
W £ C°°{U) such that a = dW . 

Notice also that, because of (jl4|) . the closeness condition da = is equivalent to a*i^ = 0. 
A straightforward consequence of the previous theorem is: 
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Corollary 3 A closed 1-form a is a solution of the Hamiltonian Hamilton-Jacobi problem if, 
and only if, d{a*H) =0. ■ 

Observe that, if a E ^^{Q) is a solution of the Hamiltonian Haniilton-Jacobi problem, as 
a*ijo = 0, then Ima is a Lagrangian submanifold of {T*Q,uj), contained in a level set of H, 
because (^HJ implies that d{a*H) = 0. In fact, dim Ima = n and, if j: Ima ^-> T*Q denotes 
the natural embedding, we have that j*u} = 0. Thus we recover some geometrical aspects of the 
classical Hamiltonian Hamilton-Jacobi theory. 

We can summarize the above results in the following: 
Proposition 7 Let a G ^^{Q) he a closed 1-form. Then, the following assertions are equivalent: 

1. a is a solution of the Hamiltonian Hamilton-Jacobi problem. 

2. d{a*H) = 0. 

3. Ima is a Lagrangian submanifold ofT*Q invariant by Zh- 

4- The integral curves of Zh with initial conditions in Ima project onto the integral curves 
of X = Tttq o Zh o Oi. 

If moreover a = dW , then these conditions can also be written as 

5. H o dW is locally constant. ■ 



Coordinate expressions Let us see how all the objects presented appear when we consider 
coordinates (g*) on Q, and the corresponding natural coordinates (g*,g*) and {q^,Pi) on its 
tangent and cotangent bundles. 

First, the coordinate expression of the Hamiltonian vector field Zh is given by 

ZH{q,p) = {q,p;dH/dp,-dH/dq). 

Consider a 1-form a G Q^{Q) and a vector field X £ X{Q). In coordinates they read 
a = ai dq^ and X = X"^ d/dq^ . Then Ta o X and Zh ° ol are vector fields along a, which read 

{TaoX){q) = (^q\a,{q);X\q),^X\q) 

(r)H BH \ 

q', aiiq); Q—ii^ "(9))' ~'^(^' "(^)) ) • 

Therefore their difference is 

{TaoX-Zno a){q) = (^q\ a^iq); X\q) - ^M,a{q)), ^^^iq) + ^iq,<q)) ] ■ 
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Its vanishing determines X = 9'H o a, that is: 

X\q) = —{q,a{q)). (15) 

With this X, the preceding difference becomes 

{TaoX-Zuo a){q) = [q\ a,(g); 0, ^{q)^{q, a(g)) + ^{q, a{q))^ . (16) 

So, the condition for a to be a solution of the generahzed Hamiltonian Hamilton-Jacobi problem 

is 

dai , , dH , , , , dH , , , , 

Now let us consider da = -rr-^dq^ A dq^ . Then i{X)da = X^ I —-4 — —4 | dq^ . On the other 

oq^ \ oq^ oq^ J 

BM f)J-f / f)J-f BM Bn \ 

hand, dH = -rr-dq' + -rr-dpi, so a*{dH) = ( -Trj{q,a{q)) + —-^{q,a{q)) j dq\ Therefore 

we have 



dq\ 

p=a{q) 



Again, with X given as 3'H o a, this expression becomes 

dq\ (17) 



.,,^, , *, ,-r-r-. fdHBai BH 

i{X)da + a*{dH) = J—J— + J— 

\BpjBq3 Bq'y p^^(g) 

Finally, if a = dW, then Oj = dW/Bq^, and the last condition in Proposition dreads 

H{q\BW/Bq') = const, 

which is the classical form of the time-independent Hamilton-Jacobi equation. 

A careful look at the local expressions (|16|) and H17|) gives an alternative understanding of 
Theorem 12 Using the vertical lift map again, we have: 

Proposition 8 Let vl he the vertical lift map of the cotangent bundle T*Q. Given a 1-form a 
on Q, and the vector field X = "JH o a on Q, we have 

vl (a, i{X)da + a* (dH)) =TaoX- Znoa. (18) 

Therefore, Ta o X — Zh o q; vanishes if, and only if, i{X)da + a*{dH) also does. ■ 

3.2 Complete solutions 

As in the Lagrangian case, we are interested in finding not only a particular solution as described 
in the preceding section, but a complete solution to the problem. In this way, we define: 

Definition 3 Consider a solution a\ depending on n additional parameters A G A, where A C 
M" is an open set, and suppose that the map <5 : Q x A ^ T*Q given by ^{q, A) = a\{q) is a local 
diffeomorphism. In this case {a\]X G A} is said to be a complete solution of the generalized 
Hamiltonian Hamilton-Jacobi problem. 
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From the definition it follows that a complete solution provides T*Q with a foliation trans- 
verse to the fibers, and that the Hamiltonian vector field Zh is tangent to the leaves. 

If {ax; A G A} is a complete solution, the integral curves of the vector fields provide all the 
integral curves of the Hamiltonian vector field Zh- This means that, if {qo,Po) G Ima, then 
there is Aq G A such that aAo('7o) = Po, and the integral curve of Xx^ through qq, lifted by ax^ 
to T*Q, gives the integral curve of Zh through {qo,po). This justifies the name of "complete 
solution" . 

Furthermore, different transversal foliations of {T*Q,uj), with Zh tangent to the leaves, are 
different ways to collect integral curves of Zh smoothly and such that they project onto Q in 
a coherent way: integral curves of Zh in Im ax project onto integral curves of the associated 
vector field Xx)- 

The relation between Zh and complete solutions is the same as in the Lagrangian case, using 
first integrals of Zh and the vector fields Xx associated to ax, for A G A. 

Finally, if the foliation is Lagrangian in {T*Q,uj), then we have a complete solution of the 
Hamiltonian Hamilton-Jacobi problem. In this case the above family of first integrals are in 
involution. 

3.3 Equivalence between the Lagrangian and Hamiltonian formulations 

This section is devoted to the equivalence between the Lagrangian and Hamiltonian Hamilton- 
Jacobi theory. We have the following: 

Theorem 3 Let {TQ,ujl,El) be a hyper-regular Lagrangian system, and {T*Q,uj,H) its as- 
sociated Hamiltonian system. Then there exists a bijection between the set of solutions of the 
(generalized) Lagrangian Hamilton-Jacobi problem and the set of solutions of the (generalized) 
Hamiltonian Hamilton-Jacobi problem. This bijection is given by composition with the Legendre 
map: X h^ a = 3^L o X . 

Proof Suppose X £ X{Q) satisfies TXoX = TloX, and T'JLoTl = Zh°'3^L. Let a = 9"LoX, 
then 

TaoX = T{JL oX) = T-JL oTXoX = T3^L oTlo X = Zh oJL o X = Zh o a , 

hence a is a solution of the Hamiltonian problem. 

Furthermore, if E^ o X = const., and a = 3'L o X, then El o ^H o a = const., that is, 
H o a = const. 

Conversely, suppose a G Q}{Q) satisfies Ta o X = Zh ° a, and F^ o "JH = T'JH o Zh, where 
3^H denotes the fiber derivative of the Hamiltonian, which satisfies 3^H = (9"L)^^, because the 
system is hyper-regular. Let X = 9'H o a, then 

TXoX = T{-JH oa)oX = T-JH oTaoX = T'JH o Zh o a = Tlo'JH o a = Tlo X , 

hence X is a solution of the corresponding Lagrangian problem. 
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In addition, if ii^ o a = const., and X = 3'H o a, then H o !JL o X = const., that is, 
El o X = const. ■ 

This result can be extended to complete solutions in a natural way. 

As is obvious, for regular but non hyper-regular Lagrangians, all this holds only in the local 
open sets where S'L is a diffeomorphism. 

3.4 Alternative Lagrangian descriptions 

Let us consider now the case of a regular system admitting alternative Lagrangian descriptions; 
that is, suppose there are regular Lagrangians functions L,L' £ C°°{TQ), L ^ L' , giving rise to 
the same dynamical vector field T^ = Tii G X{TQ) solution of the equations 

i{T)uL = dEi , i{T)uJL' = dEy . 
So, we have the same dynamics, but two different symplectic structures. 

• If X is a solution of the generalized Hamilton- J acobi problem for one of the Lagrangians, 
then it is also a solution of the generalized problem for the second Lagrangian. A sim- 
ilar result does not hold for the non generalized problem, that is, a solution X of the 
Hamilton-Jacobi problem for one of the Lagrangians will not be (in general) a solution of 
the Hamilton-Jacobi problem for the other Lagrangian, that is, X*loi = 0, but X*ijjLi ^ 0. 



• 



It is natural to compare with the situation in the Hamiltonian formalism, that is, in 
T*Q. Instead of ujl and ujl' , we have a symplectic structure u; G Q'^{T*Q), but different 
Hamiltonian functions H,H' G C°°{T*Q). Thus the same solution X of the Lagrangian 
Hamilton-Jacobi problem leads to two solutions a and a' of the Hamiltonian Hamilton- 
Jacobi problem corresponding to the Hamiltonians H and H' respectively. Nevertheless, 
notice that 



dq' dH , , dH' , ^ 



. aw' 

- dq 



where a = dW and a' = dW' locally. These are the equations for the integral curves of 
X. However, the corresponding dynamical vector fields on T*Q related to X are different. 
In other words, the difference between Zh ° a and Zh' o a' is a vertical vector field, which 
in general does not vanish. 

• The case of dynamical systems described by two alternative equivalent Lagrangians mo- 
tivates the use of two different symplectic structures in T*Q in the following way. Let 
L and L' be equivalent hyper-regular Lagrangians, and u>l and lol/ their corresponding 
Lagrangian 2-forms. If ujq £ Q,'^(T*Q) is the canonical 2-form in T*Q, we have that 
3^L*u>o = u)L- Then, let toi G Q'^{T*Q) be another symplectic structure in T*Q such that 
3'L*L0i = ujl'. Hence, 9"L' o 3^L^^ is a base-preserving transformation from {T*Q,ujo) to 
{T*Q,uJi). These transformations have been called fouling transformations j28j . 

As an example, it is not difficult to show this construction for the two-dimensional harmonic 
oscillator described both by L = ^{{v'^)'^ + (v^)^ - (g^)^ - (g^)^), and L' = v-^v'^ - q^q^ (see 
example 2 in Section [2. 2 j) . 
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These considerations show that, in order to incorporate alternative Lagrangian or Hamilto- 
nian descriptions in the Hamilton-Jacobi setting, we must introduce generahzed solutions. 

4 The Hamilton— Jacobi problem for unconstrained singular La- 
grangian systems 

In this section we are going to show how our procedure can be extended to singular Lagrangians 
without secondary constraints. 

4.1 Lagrangian formulation 

Now we consider a singular Lagrangian L € C°°{TQ). We recall that the Euler-Lagrange 
equation for a SODE T is the equation 

i{T)uJL = cIEl. (19) 

A curve 7 : M — > Q is a solution of the Euler-Lagrange equation if ^ = 7 satisfies 

i{i)uJL = (IEl o e. (20) 

We shall only consider the case of singular Lagrangians for which the following assumption 
holds: 

Assumption 1 The Lagrangian dynamical equation iiy\) has a SODE solution T G X{TQ) ev- 
erywhere defined in TQ and the rank of TUL is constant. 

The constancy of the rank of TS'L is equivalent to saying that ujl has also constant rank, 
hence lol is a presymplectic form. 

Under this assumption, the set of SODE solution vector fields is the set of sections of an affine 
bundle A -^ TQ, modeled on the vector bundle Kerr3"L -^ TQ. More precisely, the fiber oi A 
at w G TQ is 

A^ = {V£ T^{TQ) I Ttq{V) = V and i{V)ujL\v = dEiU } • 

A curve 7 : M — > Q is a solution of the Euler-Lagrange equations if, and only, the curve ^^ = 7 
satisfies £,{t) G ■Ag(i), for every t G M. Observe that if a vector V G Ty{TQ) satisfies Trqiy) = v 
then the linear 1-form i{V)uJL\v — dEilv is semibasic. 

The generalized Lagrangian Hamilton-Jacobi problem for these kinds of Lagrangians can be 
stated as follows: 



Generalized Lagrangian Hamilton Jacobi problem for unconstrained singular La- 
grangians To find a vector field X : Q ^ TQ such that, if ^■. M ^ Q is an integral curve of 
X, then ^ = X o 7: I ^ TQ is a solution of the Euler-Lagrange equation h2U\) . 

In a similar way to the first part of Section [2. 11 we have the following result: 
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Theorem 4 The following conditions for a vector field X G X{Q) are equivalent: 

1. X is a solution of the generalized Lagrangian Hamilton- Jacobi problem. 

2. X satisfies the condition lm.{TX o X) C A|iniX- 

3. X satisfies the equation i{X){X*uJL) = d{X*EL). 

4- For every v € ImX there exists w G A^j such that w is tangent to lva.X. 

5. The submanifold Im X is such that for every initial condition v G Ini X there is a solution 
of the Euler-Lagrange equations which is entirely contained in the submanifold ImX. 

Proof 

[(1) =^ (2)] Let q £ Q he arbitrary, and consider the integral curve 7 of X such that 7(0) = q. 
Denote by ^ = 7 the tangent hft of 7. Since condition (1) is satisfied and 7 is an integral curve 
of X, we have that ^(t) € -^({t)- But since .^ = X o 7 we have ^ = TX 07 = TX 0X07, which 
at t = gives TX{X{q)) = ^(0) G yi^(o) = >lx(,)- 

[(2) =^ (3)] Assume that Im (TX o X) C ^|imX) so that for every q £ Q we have that 
TqX{X{q)) G ■Ax(g)- Then, from the definition of A we have that ujL{TqX{X{q)),W) = 
{dEL\x(q),W) for every W £ Tx(q){TQ). In particular, if we take W = TgX{w) for arbi- 
trary w £ TqQ, we have that u)L(TqX{X{q)),TgX{w)) = {dEL\x(q),TgX{w)), or in other words 
{X*u;L)g{X{q),w) = {X*{dEi,)\g,w). Since this equality holds for every w G TgQ and every 
q £ Q,we deduce that i(X)(X*(jL) = X*dEL- 

[(3) =^ (4)] Let X be a vector field satisfying condition (3) and t; G ImX, so that v = X{q) for q = 
tq{v). The vector w = TX{v) satisfies the required properties. Indeed, on one hand it is clear 
that w is tangent to the image of X, and on the other we have that Ttq{w) = Ttq{TX{v)) = 
T{tq o X){v) = V, so that we have to prove that that linear 1-form i{w)uJL\v — dE^l^ vanishes. 
Since such 1-form is semibasic, we just need to prove that it vanishes when applied to elements 
of the form TX[u) for u £ TgQ: 

{i{w)ujL\v - dELU){TX{u)) = uJL\x(q){TX{v),TX{u)) - dEL\x(q){TX{u)) 

= {X*iOL)q{v,u)-d{X*EL)g{u) 
= {X*UL)g{X{q),u) - d{X*EL)g{u), 

which vanishes in view of the condition i{X)X*ujL — d{X*EL) = 0. 

[(4) =^ (1)] Assume that for every element v £ ImX there exists w £ A^, which is tangent 
to ImX. In other words, for every q £ Q (and hence v = X{q)) there exists w £ Ax(g) such 
that w = TgX{z) for some z £ TgQ. But the first condition for the element w to be in A 
is Ttq{w) = ttq{w), which for w = TgX{z) is just z = X{q). Therefore, the vector w is 
w = TgX{X{q)) and it is Ax(g)- Since this is true for every q £ Q, we have proved that 
Im(TX o X) C ^limX) which was shown to be equivalent to condition (1). 

Finally, (4) and (5) are clearly equivalent, and both are equivalent to the integrability of the 
restriction of yi to ImX (see the remark below this proof). ■ 
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Remark: Let us recall a few facts from the theory of implicit differential systems, in par- 
ticular, when the implicit system is just an affine subbundle of the tangent bundle. Let A ^ M 
be an affine subbundle of TM and let A^ be a submanifold of M. Consider the restriction A\j\f 
of the subbundle yi to A^. The following properties are equivalent: 

1. The restriction of ^l to A^ satisfies the integr ability condition for implicit differential equa- 
tions pnrr^ . 

2. For every initial condition m, £ N there exists a curve solution of the system which is 
entirely contained in N. 

3. For every m £ M there exists w € Am such that w is tangent to A''. 

Roughly speaking, the proofs of these facts are as follows: (1) and (2) are equivalent by definition 
of an integrable implicit differential system. [(2)=^(3)] is obvious: given m, £ M take the solution 
7(t) passing through m and contained in M, and then w = 7(0) is tangent to M. [(3)^(2)] 
Take a local section of ^l n TN, and an integral curve of such a section is a curve contained in 

M. 

As in the regular case, we can state the following particular problem: 

Lagrangian Hamilton Jacobi problem for unconstrained singular Lagrangians To 

find solutions X to the generalized Lagrangian Hamilton- J acohi problem for unconstrained sin- 
gular Lagrangians satisfying X*ujl = 0. 

The main results for this situation are summarized in the following: 

Proposition 9 The following assertions for a vector field X G X{Q) are equivalent: 

1. X is a solution of the Lagrangian Hamilton- J acobi problem. 

2. ImX is an isotropic submanifold of (TQ,ujl) and lui{TX o X) C ^|imX- 

3. d{X*9L) = and d{X*EL) = 0. 

4. ImX is an isotropic submanifold of {TQ, ojl) and for every v £ ImX there exists w £ Ay 
such that w is tangent to ImX. 

5. ImX is an isotropic submanifold of (TQ,u!l), and for every initial condition in IvnX there 
exists a solution of the Euler-Lagrange equations entirely contained in ImX. 

Proof They are consequences of the last theorem, taking into account that Im X is isotropic 
if, and only if, X*u)l = 0, and this is equivalent to d{X*9L) = 0. ■ 
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4.2 Hamiltonian formulation 

When the Lagrangian is singular, in general, there is no satisfactory Hamiltonian formalism 
unless certain regularity conditions hold. We will assume in what follows that: 

Assumption 2 The Lagrangian L is almost-regular, that is: P = 3^L{TQ) is a closed sub- 
manifold ofT*Q, 7L is a submersion onto its image P, and the fibers 3'L~^{3'L{p)), for every 
p G TQ, are connected submanifolds ofTQ. 

The natural embedding of P into T*Q will be denoted jq: P ^^ T*Q. Denote by 3'L^ the 
map 3'L^ : TQ -^ P defined by the relation jo ° 3'L^ = 3'L. 

For an almost-regular Lagrangian system {TQ,L) there exists a Hamiltonian formalism. 
The associated Hamiltonian system is {P,u!o,Hq), where wq = JqUJ is a presymplectic form, and 
^0 S C°°{P) is the Hamiltonian function, defined by the equation 3'L^*Hq = E^. 

The Hamilton equation is the presymplectic equation 

i{Z)uJo = dHo, (21) 

for a vector field Z £ X{P). Under our assumptions this equation has solution everywhere in 
P, although it is not unique [HI El HHI- The set of solutions is the set of sections of an afRne 
subbundle 23 ^ P of T{T*Q), modeled on the vector subbundle Ker(wo) -^ P. The fiber over a 
point a £ P is 

'Ba = {Ve T^{T*Q) I i(y)a;o|« = dHo J. 

A curve fj,: M — > T*(5 is a solution of the Hamilton equations if it satisfies 

i(/i) loq = dHo o fj,. (22) 

Hence, the curve // is a solution of the Hamilton equation if, and only if, /t(t) € '^^(t) ■ 

Bearing in mind the above comments and the results for the Lagrangian and the Hamiltonian 
regular cases, the generalized version of the Hamiltonian Hamilton-Jacobi problem for these 
kinds of singular systems can be stated in the following way, which is not exactly as in the 
regular case: 

Generalized Hamiltonian Hamilton Jacobi problem for unconstrained singular La- 
grangians To find vector fields X : Q ^ TQ such that, if j: M ^ Q is an integral curve of X 
then fi = 3'L^ 0X07 is a curve solution of the Hamilton equation \2i3^) . 

Observe that a = 3^L^ o X: Q ^ P is a section of the projection vtq = vtq o jq : P ^ Q. We 
will say that a is the 1-form associated with the particular chosen solution X . 

In this way, all the definitions, results and comments stated in Section|21hold for the manifold 
P instead oi T*Q. In particular: 

Theorem 5 The following conditions for a vector field X G 3C{Q) are equivalent 
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1. X is a solution of the generalized Hamiltonian Hamilton- J acobi problem for the uncon- 
strained singular Lagrangian L, with associated 1-form a. 

2. X satisfies the condition Im(TQ o X) C 'B|imQ, 

3. X satisfies the equation i{X){a*ijjQ) = d{a*HQ). 

4- For every A G Ima there exists w £ lix such that w is tangent to Iina. 

5. The submanifold Ima is such that, for every initial condition in Ima there is a curve 
solution of the Hamilton equations which is entirely contained in the submanifold Ima. ■ 

As above, we can state the particular case: 

Hamiltonian Hamilton Jacobi problem for unconstrained singular Lagrangians To 

find solutions a of the generalized singular Hamiltonian Hamilton- J acobi problem for uncon- 
strained singular Lagrangians satisfying a*u)Q = 0. 

And we obtain: 

Proposition 10 The following assertions for a 1-form a: Q ^ P C T*Q are equivalent: 

1. a is a solution of the Hamiltonian Hamilton- J acobi problem for the unconstrained singular 
Lagrangian L. 

2. Ima is an isotropic submanifold of P,uJo) andIm.{Ta o X) C Sjima- 

3. d{jo o a) = and d{a*H) = 0. 

4- Ima is an isotropic submanifold of {P,ujq) and, for every A S Ima, there exists w G 'Bx 
such that w is tangent to Ima. 

5. Ima is an isotropic submanifold of {P,loq), and for every initial condition in Ima there 
exists a curve solution of the Hamilton equations entirely contained in Ima. ■ 

As a final remark, the equivalence between the Lagrangian and Hamiltonian Hamilton— 
Jacobi problem in the unconstrained singular case is straightforward, taking into account how 
the problem has been stated in the Hamiltonian formalism. 

5 The Hamilton— Jacobi problem for time-dependent regular sys- 
tems 

5.1 The extended homogeneous Lagrangian formalism 

The geometric formalism for non-autonomous Lagrangian and Hamiltonian systems exhibits 
some differences with respect to the autonomous formalism (see |14[ I15| for the details) . In the 
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non-relativistic Lagrangian formalism of time-dependent dynamical systems, the configuration 
space is a bundle vr: S — > M, known as the configuration bundle, and the velocity-phase space 
is the first-order jet bundle vr^ : J^vr -^ E. Given fibered coordinates (t, g*) on E, we get fibered 
coordinates {t,q^,v'^) on J^vr. 

A non-autonomous Lagrangian is a function L G C°°(J^7r). In this case, the associated 
Cartan 1-form @l is the 1-form on J^vr whose coordinate expression is 

Q^= (dq' -v'dt)+Ldt. 

The Lagrangian is regular if the dimension of the kernel of the Cartan 2-form 0^ = —dQi is 1- 
In this case a unique vector field T, the dynamical vector field, is determined by the dynamical 
equation i{T)QL = 0, together with the normalization condition i(T)dt = 1, which ensures that 
integral curves of T are parametrized by the time coordinate t. 

It follows from the above description that our Hamilton-Jacobi theory does not apply directly 
to time-dependent systems. A way to solve this problem is to describe non-autonomous systems 
by the so-called homogeneous formalism (see P2_. for a friendly introduction see Section 2.3.1 
in JZl)) as we are about to explain. Instead of the first jet bundle J^vr, we consider the tangent 
bundle te '■ TE -^ E, which is called the extended Lagrangian phase space, and we will define 
a new Lagrangian in this extended space whose solutions are related to the solutions of the 
original system. Natural coordinates on TE will be denoted by {x^ , x^ , w^ , w'^) . 

The manifold J^vr can be canonically embedded into TE by means of the map i : J^vr -^ TE 
given by i(i/(,o") = <J(to). In fact, the image of i is included into the open submanifold TE C TE 
of vectors which are not vertical over M. Conversely, we can define a map p: TE -^ J^ir, 
which is a left inverse of i, defined as follows: for w = 7(0) G TE, we consider the function 
(f = IT o ^■. M^M; this function is locally invertible in a neighborhood of s = since 93(0) 7^ 0, 
thus we can consider £7 = 70 ip~^, which is a local section of vr. The 1-jet of a at the point 
to = ^(0) is well defined (it does not depend on the choice of the curve 7 that represents w), 
and we define p{w) = jIq'^- In coordinates, the expression of i and p are 

i(t, q\v') = (t, q\ 1, v') and plx"^ ,x\w^ ,w') = f x°, x\ ^ 

From the above expression it is clear that p o z = id (but i o p 7^ id), and it is easy to see that 
the kernel of Tp is generated by the Liouville vector field A on TE restricted to TE. 

Next we define a Lagrangian function L G C°°{TE) in such a way that the action defined by 
a curve in the jet formalism and the action defined by the corresponding curve in the extended 
formalism coincide, that is, with the same notation as above we look for a function L such that 



/ {j^a)*Ldt= / ^*Lds 

Jtn Jsn 



'to -J So 

under the change of variable t = '^{s). It follows that the Lagrangian L is defined by 

L(7(0)) = LUlaMO). 
In other words L = w^{p*L), which in coordinates reads 

L{x^,x\w^,w') = L{x^,x\w^/w^)w^ 
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which is homogeneous of degree 1. 
Proposition 11 The following relations hold: 

1. ei=p*QL andeL = i*9i. 

2. uj^ = p*VLl and Ql = ^*^i- 

3. E^ = 0. 

Proof A simple calculation in coordinates shows that 



dx'' + ^—dx' = -p*ELdx'' +p* {^—]dx'=p*{ -Eidt + T—dq' = p*Ql- 



and 

El = (Aw°)p*L - w'^p*L = w°p*L - w'^A{p*L) - w°p*L = 0. 

The other properties follow easily form the first one by using p o i = id. ■ 

When the Lagrangian L is regular, it follows from the preceding proposition that the kernel 
of uj^ is 2-dimensional and that A is in the kernel. Thus L is a singular Lagrangian, but it 
is easy to see that its dynamical equation has solutions defined everywhere and furthermore, 
among them, there are SODE solutions, since L is a type II-Lagrangian (see jS] for the details). 

d ■ d ■ ■ d 

Locally, if F = — — h v^-—^ + P{t.,q' -.v^)-—^ is the solution of the dynamics in the time- 
ot oq^ ov^ 

dependent formalism, then the solution of the dynamics in the homogeneous extended formalism 
is 

f = -°^ + -'^ + {wrf{x',x^,w^/w')^ + XA, 

for A an arbitrary function on TE. 

Once we have transformed our time-dependent problem into an autonomous one, and taking 
into account that the Lagrangian is singular but does not generate constraints, we can apply 
the theory that we have developed in the previous sections. We look for a vector field Y on E 
such that its integral curves are also integral curves of the dynamical vector fields T. Since we 
are interested in integral curves parametrized by time, we must chose such vector field Y in the 
image of the map i, that is, we will take a jet field X : E ^ J^n and the vector field Y = io X. 
Then we have that 

Y*Q^ = Y*p*eL = X*i*p*eL = X*Ql, 
so that the Hamilton-Jacobi equation amounts to d{X*QL) = 0. Locally, the form X*Ql will 
be exact, X*Ql = dS, i.e. 

dL , A dS , dS 



Thus the Hamilton-Jacobi equation reads in coordinates 

— = -EL{t,q\X') 
^ = ^U,q\X'), 
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which are the expected expressions of the Hamilton-Jacobi theory for time-dependent La- 
grangian systems. 

5.2 The Hamiltonian formalism 

In time-dependent non-relativistic Hamiltonian Mechanics, the Hamiltonian is not a function 
but a section /i of a certain bundle. Given a bundle it: E -^ R we consider the affine-dual 
bundle Aff(J^7r,M), which is canonically isomorphic to T*E, and the vector bundle u: J^*tt = 
Ver(7r)* -^ E dual to the vertical bundle. We have an affine bundle fibration /u: T*E — > J^*7r 
and a Hamiltonian is a section h of the projection ^. 

Given a Hamiltonian section h: J^*ir -^ T*E, the pullback by h of the canonical symplectic 
form u = —dO on T*E defines a 2-form n^ = h*ijj on J^*7r. The associated Hamiltonian vector 
fields are the solutions T^ to the equations 

i{Th)9.h = and i{Th)dt = 1. (23) 

It is clear that Vt^ = —dQh where Q^ = h*6. 

The relation with the Lagrangian formalism is as follows (see for details). From the 
Lagrangian L we can define two maps, usually called the Legendre transformation 9"L: J^vr — > 
J^*7r and the extended Legendre transformations 9"/, : J^vr -^ T*E, related by ^ o 9"^ = 9"L. 
When the Lagrangian is hyper-regular we have that 3'L is invertible and a unique section h of 
// is determined by the equation 7l = ho S^L. When L is regular we must restrict the study to 
the image of 3^L. For simplicity, we will assume that the Lagrangian L is hyper-regular. 

Let us consider the homogeneous Lagrangian L E C°°{TE) and the Legendre transformation 
9^j^ : TE -^ T*E defined by L. Then the relation between the Legendre transformation 9"j^ and 
9^L is given by 3"^ = z o 3"^,. In coordinates {x^ , x^ , w^ , w^) in TE and {x^,x^,u,pi) on T*E, the 
expression of the Legendre transformation is 

3^^ix^,x\w^,w') = [x\x\-p*ELy (^1^)) , 
the composition 9"^ = 3"^ o i : J^vr -^ T*E is given by 

/ PiT 

^L{t,x\v')=(t,x\-EL,^^ 

and composing with the projection jjl: T*E -^ J^*tt we get the map 3^L: J^ir -^ J^*7r, which in 
coordinates reads 

7L{t,x\v')= t,x' 



Since we are assuming that the Lagrangian L is hyper-regular, it follows that the Lagrangian 
L is almost-regular, and we can construct the Hamiltonian formulation. The kernel of T3^^ is 
spanned by the Liouville vector field A on TE^ and moreover we have Jf^{\w) = Jf^{w) for 
every A 7^ 0, so that the image of 3"^^ coincides with the image of Jl- Since JL is invertible, we 
can identify the image of 3"^ with J^*tt, or better, with the image of J^*tt by a unique section 
h: J^*7r -^ T*E of /x given explicitly hy h = 3^l o 3^L~^. 
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Thus, with the same notation as in the general case, we have that P = J^*it and jo = h. 
If we denote by w = —dO the canonical symplectic form on T*E, then the 2-form wq = Jq^ is 
ri/j = h*uj, that is the differential of Q^ = h*6. Following our general theory for unconstrained 
singular systems, we must look for a section a such that a*Qh is locally an exact form dS. In 
coordinates h{t,x'^,pi) = {t, x^ , — H {t, x^ , pi) , pi) and hence 

dS dS 

a*Qh = aidq^ — (H o a) dt = -^r-dt + -—^dx^, 

at ax* 



from where we get 
or equivalently 



ai = —^ and — + ii(t,x ,«,) = 0, 



dS ^(, i dS\ ^ 



which is the classical time dependent Hamilton-Jacobi equation. 

6 Distance on a Riemann manifold: the free relativistic particle 

6.1 General features 

We consider a Riemannian or semi-Riemannian manifold {Q,g) and the Lagrangian L{v) = 
y^g{v,v). In the semi-Riemannian case we restrict v to be time-like, i.e., g{v,v) > 0. In 
particular, if g is the Lorentz metric, this Lagrangian models a free relativistic particle on the 
manifold Q. 

Lagrangian dynamics The Lagrangian L is singular. In fact, it is homogeneous of degree 
one, hence the energy function vanishes identically El = 0. The Cartan 1-form is given by 



for all U € T{TQ), where v = ttq{U) and w = T^tq{U). The kernel of the Cartan 2-form lol is 
generated by the geodesic spray F and the Liouville vector field A. There exists underdetermined 
global second-order dynamics given by F + AA, for any function A G C°°{TQ). See [HI El- 
Hamilton— Jacobi equation Let X be a nowhere vanishing vector field on Q, and everywhere 
time-like in the semi-Riemannian case. From the expression of 9l above we immediately have 
that 

x*eL= , ^ =x\ 

where we have denoted by X the 1-form on Q such that (X , y) = giX, Y) for all vector fields 
Y on Q. If we define X as the unitary vector field in the direction of X, that is 

X= , ^ X, 

V9{X,X) 
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then we have that X*eL = XK 

Since the energy function vanishes identically, the Hamilton-Jacobi equation reduces to 
d{X*OL) = 0. Let us find an alternative expression for this condition in terms of the Levi-Civita 
connection associated with the metric. 

Proposition 12 A time-like vector field X € X{Q) is a solution of the Hamilton- J acohi equa- 
tion of the Lagrangian L{v) = y^g(v,v) if, and only if, 

1. the distribution X is integrable, and 

2. VxX = XX for some function A G C°^{Q). 

Proof If X is a solution of the Hamilton-Jacobi equation, i.e. d{X*9L) = 0, we have that 
dX^ = 0, hence X-^ = {X^)° is an integrable distribution. 

Observe that, for every vector field Z G 3£((5), using the Levi-Civita connection associated 
with the metric g, since this connection is torsion-free, the exterior differential can be calculated 
by skew-symmetrization of the covariant differential, and thus 

dX^ (y, Z) = VyX^ {Z) - VzX^ (y ) = ^(Vy X, Z) - g(VzX, Y) . (24) 

Then we have 

= {dX^){X, Z) = g(Vj^X, Z) - g{VzX, X) 

= g{Vj,X, Z) - ]^Vz{g{X, X)) = g{V^X, Z) 

then V^X = 0, and hence VxX = \X with A = Vx(ln yj g{X,X)). This proves the direct 
statement. 

Conversely, assume that X satisfies conditions (1) and (2). We will prove that dX = 0, so 
that d{X*9L) = 0. Taking the derivative Vx of g{X, X) we find that the function A is given by 
the relation \7x\/g{X, X) = Xy^g{X, X), from where Vj~-X = follows. 

On the other hand, if the distribution X is integrable, there exists locally a nowhere van- 
ishing function if such that d{ipX ) = 0. First we will prove that d(p = {'Vj^(p)X. Indeed, for 
every Z G X((5), 

= d{^X^){X, Z) = V^{^g{X, Z)) - Vz{vg{X, X)) - ^g{X, [X, Z]) 
= {Vj,ip)g{X,Z)-Vz^ 

where we have used that g{X,X) = 1 and 

g{X, [X, Z]) = g{X, V^Z) - g{X, V zX) = g{X, V^Z) - ^Vz{g{X, X)) = g{X, V^Z) 

Therefore Vzf = (Vj^(/9)(7(X, Z), for every Z G X((5), which proves that dip = {'Vj^ip)X. But 
we have 

= d{ipX^) = dip A X^ + if dX^ = ipdX\ 

that is dX'' = 0, and hence d{X*eL) =0. ■ 
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6.2 Alternative Lagrangian description 

It is well known [3] that L{v) = 29{v, v) provides a Lagrangian description of a free motion in a 
Riemannian manifold, that is, the geodetic spray T, which is, moreover, regular. It is interesting 
to compare the solutions of the Hamilton-Jacobi equation for this Lagrangian with the above 
one. 

The Lagrangian L is homogeneous of degree two. Therefore, we have that El = L = ^ffC^j '^)- 
The Cartan 1-form is given by 6l(W) = g{v,w) for W £ Ty{TQ), where w = TyTQ{W), so that 
X*9l = X' . Then the Hamilton-J acobi equation is 

dX^ = and g{X, X) = c, 

for some constant c > 0. For c = we have the trivial solution X = 0, and for c > we can 
rescale X to X/y/c, so that we can consider only the case c = 1, which means that we can 
restrict our study to the case that X is a unit vector field. 

Proposition 13 A unit vector field X is a solution of the Hamilton- J acohi equation for the 
Lagrangian L{v) = 29iv,v) if, and only if, X is integrable and VxX = 0. 

Proof By ^, we have dX\Y, Z) = g{VYX,Z) - g{VzX,Y), so that the Hamilton-Jacobi 
equation is equivalent to 

g{VYX,Z)=g{VzX,Y). 

If we take Y = Z = X the condition is identically satisfied. If we take Y = X and Z G X , 
we have g^VxX, Z) = 0, thus the vector field X satisfies that VxX = XX for some function 
A G C°°{Q). But from the normalization condition g{X, X) = 1 we have that A = gCVxX, X) = 
^Vx[9{X, X)] = 0, so that VxX = 0. Finally, for Y,Z eX^ we have giVyX, Z) = gi^ zX, Y), 
which, as above, is equivalent to g{X^ \Y, Z]) = 0. ■ 



Remarks 

1. The condition V xX = gives the generalized solution of the Hamilton-Jacobi problem, 
and together with the integrability of X^ give the classical Hamilton-Jacobi solution. 

2. Recall that a vector field X satisfying that V xX = is called a geodetic vector field, 
and its integral curves are geodesies parametrized by arc length. If we reparametrize the 
curves we have the vector field X = fX for some function / nowhere vanishing. Therefore 
V^X = f{Xf)X, so that Vx^ = ^^ with A = f{Xf). Notice the relation between the 
unit length parametrization in the regular case with the projective theory in the singular 
case. In the regular case, the vector field X must be unitary in order to have integral 
curves parametrized by arc-length. 

3. The interpretation of the above results is (in both cases) as follows: the vector X points 
in the direction of propagation of the rays, and the orthogonal distribution to X is the 
tangent to the wavefront. Wavefronts are manifolds, so the orthogonal distribution to X is 
integrable. Furthermore the rays are the geodesies of the metric, and therefore the vector 
field X must be a geodetic vector field. 
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7 Free motion on Lie groups, rigid bodies and the electron 
monopole system 

In this section we wish to show that the notion of generahzed solution is the only one available 
in generic situations, because solutions in terms of characteristic functions are not available 
globally either for topological reasons or because of invariance requirements. We are going to 
present a simplified approach to dynamics on Lie groups (see fj)' since we wish to isolate the 
main conceptual aspects of the Hamilton-Jacobi problem on spaces with nontrivial topology, 
although parallelizable. 

7.1 Free motion on Lie groups 

By free motions on a Lie group G we mean motions associated with equations of motion analo- 
gous to the equation dP'x/dt^ = 0, which are written in some affine space. Thus for simplicity we 
consider our group realized as a group of matrices g G GL{n,'E.). The equations of free motion 
will be written as 

These differential equations admit a Lagrangian description in terms of a Lagrangian function 

L{g,g) = \TT[{g-'gf] . 
The geometrical objects associated with L are simply written 

eL{g,g) = Tr [{g~^g) (g-Ug)] , ujl = -dOi, El = L . 

We will show that every left-invariant vector field X provides us with a solution of the 
generalized Hamilton-Jacobi problem. So, let us consider X G X{-jLG. Denote by ^ the value of 
X at the identity e in the group G, that is ^ = X{e) G q. In this way, we have that X{g) = gS,, 
or g-'X{g) = ^ 

On the one hand, it is clear that the pullback of the energy is constant: 
{X*EL){g) = {X*L){g) = ^Tr[{g-'X{g)f] = ^TVl^^). 

On the other hand, the pullback of the symplectic form does not vanish. Indeed, we calculate 
{X*6l,Y) for a vector field Y £ X(G), which we may take to be left-invariant, Y{g) = g(^, for 
some C ^ Q- 



{X*9L,Y){g) = ^L{X{g)+tY{g)) = 11 Tr[(e + tC)'] = Tr(eC). 



d 
di' 



t=o 2 dt 



t=o 



Thus the differential evaluated on two left-invariant vector fields Yi , 12 , Yi {g) = g(i and I2 id) 
9C2, is 

d{X*9L)ig)iYi,Y2) =Yi{g){Tr{CC2)) -Y2{g)iTtiCCi)) -TTiC9-'[Yi,Y2]{g))) 
= 0-0-Tr(C[Ci,C2]), 
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so that {X*lol){Yi,Y2) = Tr(e [(1,(2]). 

It follows that X is not a solution of the Hamilton-Jacobi problem in the standard sense 
(except for Abelian groups). Nevertheless X is a solution of the generalized Hamilton-Jacobi 
problem, since ixiX*ujL) = 0: 

ix{X*u;L)iY) = Tr(e [t C]) = TrifC " CCO = 0, 

where use has been made of the properties of the trace. 

Therefore, the left trivialization provides a diffeomorphism <I> : G x g -^ TG, given by 
^{9^0 = ^e-^g(C)) such that for every S, £ Q we get a solution of the generalized problem. 
In other words we have a complete solution of the generalized Hamilton-Jacobi problem, where 
the parameter space A is the Lie algebra, A = g. 

The associated constant of the motion F = pr20$^^ : TG ^ g is explicitly given by F{g, g) = 
g^^g, or in other words F{g, g^) = ^. We can identify q with q* via the trace operation, that is 
we identify z^ G g* with .^ G g if {i^, (") = Tr(,^C) for every C ^ S- Under this identification, we get 
a map fj, : TG -^ g* given by 

{f,{g,g),C) = Tv{g-'gC), 

which is the momentum map for the left action of G on TG. 

Notice that we have exploited the left-invariance of the Lagrangian function. Furthermore 
L is also right invariant, since we can write it in the form L{g,g) = ^ TT[{gg^^)'^]. Therefore we 
can also define a second foliation by taking the right-invariant vector fields. 

Finally, it should be remarked that whenever G is a compact Lie group we cannot have 
functions on G whose differentials are never vanishing, therefore any invariant foliation (solutions 
of the generalized Hamilton-Jacobi problem) could never be associated with some dW, for some 
function W : G x A ^ R. 

7.2 Rigid bodies 

Consider a (generalized) rigid body defined on a configuration Lie group G with symmetric 
inertia tensor /: g -^ g*. We will analyze it in the Hamiltonian formalism on the cotangent 
bundle T*G. The Hamiltonian function H G G°°{T*G) is 

^(A,) = ^{T:L,{X,),r'T:Lg{X,)), 

where ( , ) is the standard pairing. We will show that every right-invariant 1-form a G i}^{G) 
is a solution of the generalized Hamilton-Jacobi problem. If ^u G g* is the value of a at the 
identity e G G, then we have that a{g) = T*i?„-i(/i). The value of a on a right-invariant vector 
field Y = TeRg{C) is constant, 

{a,Y){g) = {T;Rg-ifi,T,Rg{C)) = {fi,C), 

and hence the differential of a over two right-invariant vector fields Yi and I2 is 

da{Yi,Y2){g) =Yi{g){a,Y2) -Y2{g){a,Yi) - {a,[Yi,Y2]){g) 

= Yi{g){fi, C2) -Y2{g){fi,Ci) - {fi,-[CiX2]) 

= 0-0 + (m,[Ci,C2]>. 
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Therefore, for every Ci; C2 G we have 

da{g){T,Rg{Ci),T,Rg{ZC2)) = {f^,[Ci,C2]). 

The fiber derivative of the Hamiltonian is given by 3^H{Xg) = T(,LgI^^ T*Lg{Xg), for every 
X-geT*G. Indeed 

{X'g,FH{Xg)) = lH{Xg+tX'^)\^^^ = {T:Lg{X'g),I-'T:Lg{Xg)) = {X'g ,T,L gl'^ T*^ L g{X g)) . 

Therefore, the vector field X G X(G) associated to a is 

X{g) = ^H{a{g)) = 3^HiT*Rg-,{fi)) = T,Lgr\Ad*g{fi)), 

The contraction of da with X is given by 

iixda){g){TRg{C)) = da{g){X{g),TRg{C)) 

= da{g){nRg{TgRg-iX{g)),TRg{C)) 

= {fi,[AdgI''Ad*gfi,C]), 

where we have used that 

T;Rg-i{X{g)) = T*Rg-iTeLgr^Ad*gH = Adgr^Ad*gfi. 



{a*H){g) = - = {T:Lg{T;Rg-.{ii)),I-^T:Lg{T;Rg-.{,,))) = -{Ad*gip),I-'Ad;{^,)). 



{dia*H)ig),T,RgiC)) = -{a*H){jit)) 



t=o 

-1 



Furthermore, the puhback of the Hamiltonian by a is 

^ = {T:Lg{T;Rg-^{i,)),r^T:Lg{T;Rg-r{i,))) = I 

To calculate its differential evaluated at TeRg{C), we consider its integral curve 7(t) = eyip{t()g 
through the point g and hence 

d_ 

= {Ad*gadlii, r^Ad*gH) 
= {adlfi, Adgl^^ Ad*gix) 
= {fi, ad(^AdgI^^ Ad* fj) 

and finally adding both terms we get 

{i{X)da + d{a*H))(TeRg{C)) = (^, [Adgr^Ad*gfi,C]) + {f^,adt;Adgr^Ad*gfi) = 0. 

Thus we have a complete solution of the generalized Hamilton-Jacobi problem, <1> : G x g* ^ 
T*G explicitly given by the inverse of the right trivialization map, ^{g,fi) = T*Rg-i{fi). The 
associated constant of the motion is F = pr2 o <|)^^ : T*G -^ g*, which is the momentum map 
F = Jl associated to the left action of G on T*G, that is, 

F{Xg)=T:RgiXg). 

As the theory predicts, if g{t) is an integral curve of the vector field X, i.e. it satisfies 
g{t) = TeLg(£^I~^{Ad*,^.^), then Q. = g^'^g is given by 10. = Ad*iJ. (with ^ G s* constant) and 
hence it satisfies the differential equation IQ = —ad^Ad*^ = —ad^{10), which is the Euler 
equation. 
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7.3 The electron monopole system 

The equations of motion for a charged particle with electric charge e moving in the external 
magnetic field of a monopole with magnetic charge g are described by the following second order 
vector field in Q = M^ - {0} 



d n j k ^ 

where e^-fc is the completely skew-symmetric Levi Civita tensor, i.e. such that ei23 = 1 and with 



Qqj J.3 *J Qyt 



eg 
n 



4:TTm 



The vector field T admits a symplectic description with the symplectic structure (see e.g 

m) . . n ■ ■ 

u = dx^ A dv^ 5- Eijk x^ dx^ A dx^ , 

2fi ■> 

and Hamiltonian function 

H = —vP rP . 

Because uj is closed but not exact, uj cannot be written as a Lagrangian 2-form uji^. It is however 
possible to write it as a Lagrangian 2-form locally by using a local Lagrangian. 

In addition to the Hamiltonian function, the dynamical system admits other constants of 
the motion associated with the rotational symmetry group; they are 

/A, 

n — fijfc X V -\ . 

They are made up of the expected components of the orbital angular momentum plus the 
"helicity term" nx^ jr. 

It is possible to find local solutions of the standard Hamilton-Jacobi equation by using con- 
stants of the motion //, 1 and /a, for instance. We may solve for the velocities, and by replacing 
them in L we find a 3-parameter family of vector fields defined on some open submanifold of 
M^ - {0}. 



It should be noticed, however, that it is not possible to find globally defined vector valued 

d 



solutions, because if we denote the sought solution by y = Y^ jfj, we would have 



Y*{dx' A dv"-) = — -ejki x^ dx^ A dx" , 
2 r'^ 

which is not possible because the left hand side is exact while the right hand side is not. 

Nevertheless, it is possible to describe the electron monopole system as a reduction of a glob- 
ally defined Lagrangian system with a singular Lagrangian but without secondary constraints. 
To this end we replace the configuration space Q = M^ — {0} ~ 5^ x M-|_ with a covering 
by replacing S^ with S^ in the product of manifolds. The new configuration space will be 

5^7(2, C) xR+. 

The covering map vr : SU{2,C) — > S'^ is given by the following construction. Let {x^,x'^,x^) 
be the coordinates in M.^ — {0} and let x^ = x^ /r S 5^, so that they satisfy x^ x^ = 1. Now we 
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describe 
we have 



in terms of the 2x2 traceless Hermitian matrices using as a basis Pauh matrices, 

„3 



M = x-a 



X 

x^ + ix'^ 



1 • 2 

X — IX 



-X 



Now we describe our covering map by introducing the following matrices to describe 

s = y^ I + iy ■ a ■ 



yU _|_ ^yO y^ _|_ ^yl 

o .1 n . Q 

—y + ly y — ly 



and setting vr : M^ — > M'^ by means of 



7r(s) = sa^s^ = X ■ a . 

This map is also known as the Kustaanheimo-Stiefel map (for a classical and quantum version 
of this map see the recent papers jUISI)- 

This relation makes sense because both sides are traceless Hermitian matrices and Pauli 
matrices are a basis for the real linear space of Hermitian matrices with zero trace. We notice 
that s represent elements of SU{2, C) when the constraint 



ssT = {{yy + (y^ + (y^ + (y^)^) / = det s I = I 



is imposed. 



To spell out the way {x^,x'^,x^) depends on (y^ , y^ , y'^ , y^) , i.e. the pull-back of coordinate 
functions from M^ — {0} to M.^ — {0}, we notice that 



.t 



yU _ ^yO _y^ _ ^y± 

y^ - ty^ y^ + ly^ 



so that 

sass 

which is given by 



y" _|_ ly-^ y^ _)_ ly^ 

9 .1 n . o 

—y + ly y — ly 



1 
-1 



n • s 2 -1 

y^ — ly^ —y — ly 

y — iy y + iy 



S (73 s^ = 
provides us with 



(y0)2 ^ (^3)2 _ (yl)2 _ (y2)2 2(y0 yl + y2 y3) _ 2i(y0 y2 _ yl ^3) 
2(y0 yl + y2 y3) + 2i(y0 y2 _ ^l ^3) .(^0)2 _ (^3)2 ^ (^1)2 ^ (y2)2 



x^ = 2(y%i+y2y3) 
x2 = 2(y%2_yiy3) 

3 _ /„,0\2 I ^,,3\2 /,,1\2 /„,2n2 



X 



{yy + (y-r - (yr - (yy 



Now we find the tangent map of the covering map 



Ttt : T{SU{2,C) x 1R+) -^ T{S^ x 



^+J, 



more explicitly 



2(y%i + u°yi + u^y^ + /u^) 
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v^ = 2(y°u° + u^y"^ - y^u^ - y'^u^) 
and the pull-back of the 2-form 

Lo = dx^ A df * T eijk x* dx^ A dx^ , 

namely, {Tti)*uj, will be exact and moreover Lagrangian 

{Ttt)*uj = UJL, 
with Lagrangian L on T{SU{2,C) x IR+) given by 

L = —T7r*{v^v^) + iriTvasS^ s. 

The fibering map vr : SU{2, C) x M+ ^ S"^ x ]R_|. is actually a principal bundle projection 
with group ^(1) given by 

C/(l) = {exp(ito-3) I t GM}, 

and acting on SU{2, C) on the right. The tangent bundle group TU{1) = U{1) (8> IR will now be 
the structure group of the tangent bundle T{SU{2,C) x M+) -^ T{S'^ x M+). 

Within the notation we have used, the left-invariant vector field along a^ generator of the 
U{1) action is 

X^ = yO A _ y3 A _ £3^ . y* A. 

Qy3 QyO *■? Qyj 



while the infinitesimal generator of TU{1) will be the tangent lift 



and the vertical lift 

By using the pull-back of constants of the motion fromT(]R'^ — {0}) to T(M^ — {0}), we will find 
generalized extended space. On the 8-dimensional carrier space they will define submanifolds of 
codimension three. If the constants of the motion used are pairwise in involution, the invariant 
submanifold will be isotropic, otherwise we will give rise to solutions of the generalized Hamilton- 
Jacobi problem. 



8 Conclusions and outlook 

In this paper we show that to deal with bi-Hamiltonian systems in the Hamilton-Jacobi setting 
it is convenient to introduce generalized solutions, i.e. invariant submanifolds (or foliations) 
with dimension equal to the dimension of the configuration space without the requirement of 
Lagrangianity. Thus the associated pde will have solutions given by vector valued functions. 
When the "Lagrangian" requirement is made, these functions will be the coefficients of an exact 
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1-form, and we recover the standard PDE for the principal function W or the characteristic func- 
tion S. The hnk via the Feynman approach to quantum mechanics between these solutions and 
the phase of the wave function seems to suggest that only invariant foliations with Lagrangian 
leaves with respect to the admissible alternative symplectic structure should be accepted. 

According to von Neumann's representation theory we would have to accept as Hilbert 
spaces the space of square integrable functions defined on some invariant Lagrangian submanifold 
(according to the chosen symplectic structure). This raises the problem of the selection of 
the appropriate "Lebesgue measure" and how to compare the descriptions on these alternative 
Hilbert spaces. These aspects will be taken up elsewhere. 

Formulation of the Hamilton-Jacobi theory on the tangent bundle in terms of the Lagrangian 
formalism, prepares us ready to consider the problem of Hamilton-Jacob theory connected with 
degenerate Lagrangians (gauge theories) in full generality. Extension of the ideas in this paper 
for Lagrangian and Hamiltonian systems on Lie algebroids |3H I25j is also worthy of study. These 
aspects should be addressed in the future. 

Acknowledgments 

Support of projects MEC (Spain) BFM-2003-02532, BFM-2002-03773, BFM-2002-03493, MTM2005- 
04947, and FPA-2003-02948 and C02-399 is acknowledged. We wish to thank Mr. Jeff Palmer 
for his assistance in preparing the English version of the manuscript. 



References 

[1] R. Abraham, J.E. Marsden, Foundations of Mechanics, 2nd edition, Benjamin-Cumming, Read- 
ing, 1978. 

[2] A. D'AvANZO, G. Marmo, "Reduction and unfolding : The Kepler problem". Int. J. Geom. 
Methods Mod. Phys. 2 (2005) 83-109. 

[3] A. D'AvANZO, G. Marmo, A. Valentino, "Reduction and unfolding for quantum systems: The 
Hydrogen atom". Int. J. Geom. Methods Mod. Phys. 2 (2005) 1043-1062. 

[4] P. G. Bergmann, Introduction to the Theory of Relativity, Prentice-Hall, New York, 1942. 

[5] F. Cantrijn, J.F. Carinena, M. Crampin, L.A. Ibort, "Reduction of degenerate Lagrangian 
systems", J. Geom. Phys. 3 (1986) 353-400. 

[6] J.F. Carinena, "Theory of singular Lagrangians", Fortschr. Phys. 38 (1990) 641-679. 

[7] J. F. Carinena, M. Crampin, L. A. Ibort, "On the multisymplectic formalism for first order 
field theories", Diff. Geom. Appl. 1 (1991) 345-374. 

[8] J. F. Carinena, C. Lopez, "The time-evolution operator for singular Lagrangians", Lett. Math. 
Phys. 14 (1987) 203-210. 

[9] J.F. Carinena, C. Lopez, N. Roman-Roy, "Geometric study of the connection between the 
Lagrangian and Hamiltonian Constraints", J. Geom. Phys. 4 (1987) 315-334. 



J.F. Carinena et al, Geometric Hamilton- J acohi theory 39 

[lo; 
[11 



[12; 
[13; 

[14 
[15 
[16 

[17; 
[is; 

[19 

[20 

[21 

[22; 
[23; 

[24 

[25 

[26; 

[27; 

[28; 

[29; 



J.F. Carinena, M.A. del Olmo, P. Winternitz, "On the relation between weak and strong 
invariance of differential equations", Lett. Math. Phys. 29 (1994) 151-63. 



C. C. Chiang, S. C. Lee, G. Marmo, "Lagrangian dynamics on higher-dimensional spaces with 
application to Kaluza-Klein theories", J. Math. Phys. 26 (1985) 1083-1092. 

M. Crampin, "On the differential geometry of the Euler-Lagrange equations and the inverse prob- 
lem of Lagrangian dynamics", J. Phys. A: Math. Gen. 14 (1981) 2567-2575. 

M. Crampin, "Tangent bundle geometry for Lagrangian dynamics", J. Phys. A: Math. Gen. 16 
(1983) 3755-3772. 

M. Crampin, G.E. Prince, G. Thompson, " A geometrical version of the Helmholtz conditions 
in time-dependent Lagrangian dynamics", J. Phys. A: Math. Gen. 17 (1984)1437- 1447. 

M. Crampin, E. Martinez, W. Sarlet, "Linear connections for systems of second-order ordinary 
differential equations", Ann. Inst. H. Poincarc 65 (1996) 223-249. 

B. DuBROViN, G. Marmo, A. Simoni, "Alternative Hamiltonian description for quantum sys- 
tems". Mod. Phys. Lett. A 5 (1990) 1229-1234. 

G. ESPOSITO, G. Marmo, G. Sudarshan, From classical to quantum mechanics, Cambridge 
University Press, 2004. 

Z. Ge, J. E. Marsden, "Lie- Poisson Hamilton-Jacobi theory and Lie-Poisson integrators" , Phys. 
Lett. A 133(3) (1988) 134-139. 

X. Gracia, "Fibre derivatives: some applications to singular lagrangians" , Rep. Math. Phys. 45 
(2000) 67-84. 

X. Gracia, J. M. Pons, "On an evolution operator connecting Lagrangian and Hamiltonian 
formalisms", Lett. Math. Phys. 17 (1989) 175-180. 

X. Gracia, J.M. Pons, "A generalized geometric framework for constrained systems" , DiiF. Geom. 
Appl. 2 (1992) 223-247. 

J. Klein, "Espaces variationnels et mecanique" Ann. Inst. Fourier (Grenoble) 12 (1962) 1-124 

J. Klein "Operateurs diffcrentielles sur les varietes presque tangentes" C. R. Acad Sci. Paris 257 
(1963) 2392-2394 

M. de Leon, J. Marin-Solano, J.C. Marrero, M.C. Munoz-Lecanda, N. Roman-Roy, 
"Singular Lagrangian Systems on Jet bundles", Fortschr. Phys. 50 (2002) 105-169. 

M. DE Leon, J. C. Marrero, E. Martinez, "Lagrangian submanifolds and dynamics on Lie 
algebroids", J. Phys. A.: Math. Gen. (Topical Review) 38 (2005) R241-R308. 

P. Libermann, cm. Marle, Symplectic geometry and analytical dynamics, D. Reidel Publishing 
Company, Dordrecht, 1987. 

G. Marmo, G. Morandi, N. Mukunda, "A geometrical approach to the Hamilton-Jacobi form 
of the dynamics and its generalizations", Riv. Nuovo Cim. 13 (1990) 1-74. 

G. Marmo, N. Mukunda, J. Samuel, "Dynamics and symmetry for constrained systems: a 
geometrical analysis", Riv. Nuovo Cim. 6 (1983) 1-52. 

G. Marmo, C. Rubano, Particle Dynamics on Fibre bundles. Bibhopohs, Napoh, 1988. 



J.F. Carinena et al, Geometric Hamilton- J acobi theory 40 



[30] G. Marmo, E.J. Saletan, A. Simoni, B. Vitale, Dynamical Systems: A Differential Geometric 
Approach to Symmetry and Reduction. Wiley, Chichester, 1985. 

[31] E. Martinez, "Lagrangian Mechanics on Lie algebroids". Acta Appl. Math. 67 (2001) 295-320. 

[32] G. Mendella, G. Marmo, W.M. Tulczyjew, "Integrability of implicit differential equations" , 
J. Phys. A: Math. Gen. 28 (1995) 149-163. 

[33] G. Morandi, C. Ferrario, G. Lo Vecchio, G. Marmo, C. Rubano, "The inverse problem in 
the calculus of variations and the geometry of the tangent bundle", Phys. Rep. 188 (1990) 147-284. 

[34] G.E. Prince, J.E. Aldridge, G.B. Byrnes, "A universal Hamilton-Jacobi equation for second- 
order ODEs", J. Phys. A: Math. Gen. 32 (1999) 827-844. 

[35] J. Samuel, "A phase space approach to genmeralized Hamilton-Jacobi theory", J. Math. Phys. 21 
(1980) 55-59. 

[36] W.M. Tulczyjew, "Les sous-varietes lagrangiennes et la dinamique hamiltonienne" , C.R. Acad. 
Sci. Paris 283A (1976) 15-18. 

[37] A.M. Vinogradov, LS. Krasil'schik. "What is Hamiltonian formahsm?", Russian Math. Sur- 
veys 30(1) (1975) 177-202. 



